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The spectrum of tight binding electrons on a square lattice with half a magnetic flux quantum 
per unit cell exhibits two Dirac points at the band center. We show that, in the presence of an 
additional uniaxial staggered potential, this pair of Dirac points may merge into a single one, with 
a topological transition towards a gapped phase. At the transition, the spectrum is linear in one 
direction and quadratic in the other one (a spectrum recently named "hybrid" or "semi-Dirac"). 
This transition is studied in the framework of a general Hamiltonian describing the merging of Dirac 
points. The possibility of creating gauge flelds for cold atoms in optical lattices may offer the flrst 
opportunity to observe this merging of Dirac points and the hybrid dispersion relation. 
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INTRODUCTION 



Condensed matter offers the possibility of manipu- 
lating the energy spectrum of electrons and modifying 
their free dispersion relation. Band theory leads to com- 
plex dispersion relations, commonly with quadratic ex- 
pansions in the vicinity of peculiar symmetry points. 
These quadratic expansions are characterized by a tensor 
of eflPective masses, possibly with positive and negative 
masses (near a saddle point). Graphene offers the excit- 
ing situation where the dispersion relation at low energy 
is linear^ similar to the spectrum of relativistic particles 
described by the Dirac equation. ^'^ More precisely, the 
spectrum has the form of two cones (the so-called "Dirac 
cones" or "Dirac points") in the vicinity of the K and 
K' points of the reciprocal space. It has been recently 
proposed and studied the existence of a hybrid spectrum, 
quadratic in one direction and linear in the other one, 
with interesting consequences for the energy spectrum of 
Landau levels in the presence of a magnetic field. ^ Such 
a spectrum may appear in a hypothetical graphenelike 
structure where one {t') of the hopping integrals {t) be- 
tween nearest neighbors is increased. ^"""^^ By doing so, 
the two Dirac points move and, for a critical value of this 
hopping integral {t' = 2t), they merge into a single one, 
with the peculiar hybrid spectrum. 

Other possible systems exhibiting a hybrid spectrum 
have been proposed including the organic conductor a- 
(BEDT-TTF)2l3,^^"^^ V02/Ti02 nanostructures,^'^ or 
cold atoms trapped in an optical honeycomb lattice. 
In the latter example, a recent extensive study has shown 
the possibility of tuning the position of the Dirac points 
by changing the intensity of the laser fields. In the con- 
text of V02/Ti02 nanostructures, the hybrid point has 
been baptized a "semi-Dirac" point, a name that we 
will use in this paper. Quite recently, one of us has pro- 
posed a general framework to study the motion and the 
merging of Dirac points in 2D crystals with time-reversal 
and inversion symmetries, within the framework of the 
following Hamiltonian^'^ 
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where, by varying the product m*A from negative to 
positive values, a topological transition separates a phase 
with two Dirac points from a gapped phase, with a semi- 
Dirac spectrum at the transition, i.e. when A = 0. This 
Hamiltonian is universal in the sense that the parameters 
do not refer to a specific crystal and can be related to 
the microscopic parameters of any 2D system. It has 
been shown that in a magnetic field t he L andau level 
spectrum evolves continuously from a \friB dependence 
to a linear (n + 1/2)5 behavior, with a [(n + l/2)Bf^^ 
dependence at the transition. ^'^'^ 

At the moment, there is no straightforward experimen- 
tal evidence for such merging of Dirac points in electronic 
systems. A very interesting alternative is the possibility 
of fabricating a "crystal" of cold atoms in an artificial 
optical lattice. Unfortunately, the atoms are neutral and 
are not coupled to a magnetic field. A fictitious mag- 
netic field can be induced by a rotation of the sample 
(then the role of the magnetic field is played by the Cori- 
olis forcc""^^), or by effective gauge potentials. Recently 
there have been proposals for the realization of the Hof- 
stadter spectrum in a square lattice of cold atoms, ^^'^^ 
using such U{1) gauge potentials, mimicking an external 
magnetic field. Non-Abelian gauge potentials have also 
been proposed. 

The so-called Hofstadter spectrum represents the en- 
ergy levels of two-dimensional electrons on a square lat- 
tice in a magnetic field, within the nearest-neighbor tight 
binding model. The spectrum depends on the value of 
the magnetic fiux (j) through an elementary plaquette of 
the lattice. For rational values of the reduced magnetic 
flux, (f> = = p/q^ 00 = h/e being the flux quantum, 

the electronic spectrum consists of q sub-bands, leading 
to the famous fractal structure of the spectrum. 

It is well-known that for a half-fiux quantum, (j) = 
00/2 = /i/2e, the spectrum exhibits two Dirac cones 
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at the band center, quite similarly to the spectrum of 
graphene. One may wonder if these two Dirac points can 
be manipulated (moved and merged) by some additional 
parameter, somehow similarly to the modification of the 
transfer integrals in the honeycomb lattice. In this pa- 
per, we show that an appropriate additional parameter 
is a staggered potential applied along one direction. The 
addition of this uniaxial staggered potential modifies the 
position of the Dirac points which may eventually merge 
into a "hybrid/semi-Dirac" point for a critical value of 
the potential. The structure of the paper is the following. 
In section II, we review the basic equations for electrons 
on a square lattice in a magnetic field and a uniaxial 
staggered potential. We discuss the general structure of 
the resulting Hofstadter spectrum. After a brief analysis 
of the low field spectrum (sec. Ill), we present our main 
results in section IV on the evolution of the spectrum for 
a fiux in the vicinity of 0o/2. After a short analysis of 
the spectrum near the extrema of the band, we study the 
spectrum near the center of the band (e = 0). For ^o/2, 
it consists of two Dirac cones which merge for a critical 
value of the staggered potential. Using the mapping to 
an effective Hamiltonian near e = 0, we describe quanti- 
tatively the evolution of the Landau levels in the vicinity 
of 00 /2. We conclude in section V. 



II. BUTTERFLY SPECTRUM WITH A 
UNIAXIAL STAGGERED POTENTIAL 



magnetic field 5, using the Landau gauge where Ay = 
Bx, the Schrodinger equation reads 

-^0m+l,n - ^0m-l,n + (-l)'^As0^,-^ (2) 

where 0m, n is the amplitude of the wave function on site 
{m^n), t is the hopping integral between nearest neigh- 
bors and (p = eBa^/h is the dimensionless flux through 
an elementary plaquette of the lattice. Consider the case 
of a commensurate flux (p = p/q where q is even. There 
are q inequivalent sites in the unit cell, which is now q 
times larger along the x direction. We introduce a cell 
index /, so that m = ql -\- where j is the position of 
the site in the unit cell: j = 1, • • • , Since there are q 
sites per unit cell, the Brillouin zone in the x direction is 
q times smaller and Bloch's theorem implies that 

with fe = {kx, ky) and we have 

V'r^=V'i^^ (4) 

The Schrodinger equation now reads 

M^'' = -ti^^^'-U^-'' (5) 
-2til)^^^ cos {kya - 2nj^) + (-l)^A,vi') 
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FIG. 1: Square lattice in a uniform magnetic field and a uni- 
axial staggered potential along the x direction. The black and 
white discs represent the on-site potentials =b As . ip is the di- 
mensionless flux per plaquette. 

We consider the problem of tight binding electrons on 
a square two-dimensional lattice. The sites are written 
as X = ma and y = na, a being the lattice spacing. In 
addition, we apply a uniaxial staggered potential of the 
form {—l)^As along the x direction, as illustrated in Fig. 
(1). The unit cell is made of two atoms having different 
site energies zbA^. In the presence of a perpendicular 



iqkxa 



(6) 



-2t'0^^^ cos {kya ■ 



2TTip) - A. 



(7) 



2^7/;^^^ cos {kya - 27rqip) + {-lyAs^p)^ 



(9) 



This is a g X g system, with q eigenvalues. If q is odd, 
sites with even and odd m are inequivalent, and the unit 
cell has a size 2qa in the x direction. Eqs. (6,7,8) are 
unchanged (except q ^ 2q since j = 1 , • • • , 2q) and 

/(2g+l) /(I) 

The flgures (2) present the evolution of the spectrum 
(energy versus magnetic flux) when the staggered poten- 
tial Ag is increased. Fig. (2. a) is the familiar Hofstadter 
spectrum obtained for A^ = 0. Fig.(2.b) shows the spec- 
trum for Ag = t. The most striking difference between 
(2. a) and (2.b) is that many gaps have been fllled. This 
is even more spectacular in Fig.(2.c), corresponding to 
the critical value A^ = 2t of the staggered potential for 
which the Dirac points in the if = 1/2 dispersion rela- 
tion merge, as we will see later. This suppression of 
many gaps is qualitatively understood by the the fact 
that, when A^ increases, there are more and more open 
classical orbits which are not quantized, as we discuss in 
subsection III A. 
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The main goal of our paper is the study of the spectrum 
near half flux quantum Lp = 1/2, near the center of the 
band. For Ag = 0, this spectrum is quite similar to the 
low field spectrum of graphene : due to the presence of 
two Dirac points, the spectrum near (p = 1/2 consists 
of a series of Landau levels varying as ^/nf^ where / = 
I (/p — 1/2 1 is the deviation from cp = 1/2 (Fig. 2. a). These 
Landau levels are doubly degenerate, due to the two- fold 
degeneracy of the Dirac spectrum. When increases, 
the degeneracy is progressively lifted (Fig.2.b) and, at the 
critical point Ag = 2 (Fig.2.c), the levels vary as [{n + 
1/2)/]^/^, as we will show in subsection IV C. Then, for 
Ag > 2, a gap opens and the Landau levels progressively 
evolve towards a linear variation as (n + 1/2)/. 

We now discuss the different parts of the spectrum, 
and their evolution when the staggered potential A^ is 
increased. 



III. LOW FIELD SPECTRUM 

A. Zero field spectrum 

In zero field, due to the staggered potential, the Hamil- 
tonian has a 2 x 2 structure 

"Hik) = -2 I ~ I (8) 

^ ^ y COS kx cos ky -\-r J ^ ^ 

where we introduce the dimensionless parameter r = 
As/{2t) = As/2, and where we choose from now t = 
a = h = e = 1. Therefore the zero field spectrum is given 
by 

e{k) = -2 cos ky ± 2 Vcos2"fc~+r2 (9) 

As shown in Fig. 3, the spectrum consists of two bands 
which overlap in energy. The lower band extends from 
^rnin ^ _2 _ 2^1 + to ef^'^^ = 2 - 2r, with two 
saddle points (which yields logarithmic singularities in 
the density of states) at energies e^^ = — 2 — 2r and 
= +2 — 2VI + r^. Symmetrically, the upper band 
extends from e^^*^ = -2 + 2r to e'//''^ = 2 + 2vTTr2, 
with two saddle points (with logarithmic singularities) at 
energies e^^ = -2 + 2^1 + and ef = 2 + 2r. The po- 
sition of these points is shown on Fig. (4) on a plot of the 
density of states. 

In each band, the saddle points separate regions with 
open orbits (for energies between the two saddle points) 
and closed orbits. When Ag = 0, all orbits are closed 
and the total bandwidth for incommensurate fiuxes is 
zero. This can be qualitatively understood on semiclassi- 
cal grounds: all closed orbits are quantized and magnetic 
breakdown induces tunneling between orbits which are 
all quantized. When A5 increases, there are open or- 
bits which are not quantized and the total bandwidth for 
incommensurate fiuxes increases like the number of open 
orbits. This qualitatively explains why the spectrum is 
more and more dense when Ag increases. 




FIG. 2: Evolution of the Hofstadter spectrum for the square 
lattice, in the presence of a uniaxial staggered on-site potential 
characterized by the parameter r = As/(2t). a) r = 0, this 
is the usual Hofstadter spectrum. In the center of the band 
(e = 0) near cp = 1/2^ the energy levels vary as \fnf , where 
f —\ip) — l/2\. b) r — 1/2, the degeneracy of the levels near e = 
0; ^ — 1/2 has been lifted, c) Spectrum at the critical point 
r = 1. The levels near e = 0, = 1/2 vary as [(n+ 1/2)/]^/^. 
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FIG. 3: Zero field spectrum e{k) for tight binding electrons on 
the square lattice with a staggered potential. The two bands 
overlap when r < 1. Left panel: r = 0.5. Right panel: when 
r = 1, the two bands no longer overlap. We have represented 
the plane e = . 



B. Low field spectrum, at the edges of the band 

At low magnetic field {(p <C 1), the energy levels near 
the extrema of the band, are grouped into linear Landau 
levels (Fig. 2). Their spectrum is obtained from a low 
energy expansion of the dispersion relation (9) near k — 
0: 



e(k) = t2(1 + \/l + r2) ± 



(10) 



leading to a cyclotron mass mo = (l+r2)i/'^/2 and there- 
fore to a set of Landau levels given by {B — 27r(/p in our 
notations) 



47r 



enM = T(2+2VTT72)± 



(rT7)vi(-+2)^+---(ii) 



which describes properly the evolution of the Landau lev- 
els at the extrema of the spectrum (Figs. 2). 
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FIG. 4: Zero field density of states for tight binding electrons 
on a square lattice with a staggered potential, corresponding 
to the dispersion relation (9), for r = (a), r = 0.5 (b) and 
r — 1 (c). We have plotted the density of states in the lower 
band (dashed line), the density of states in the upper band 
(dotted line) and the total density of states (full line). 



C. Low field spectrum, center of the band 

Figure (2. b) shows that, in the middle of the band and 
for a finite A^, two set of Landau levels intermix. This 
is emphasized in the enlargement displayed in Fig. (5). 
This structure is easily understood by looking at the zero 
field spectrum, displayed on Fig. (3). When < r < 1, 
the spectrum is formed of two overlapping bands. The 
maximum of the lower band is located at energy e = 
2(1 — r) and the minimum of the upper band has energy 
e = — 2(1 — r). An expansion around these energies gives 

e(g±) = T2(l-r)±(fc^ + M) + ... (12) 



where gr^ are deviations respectively to the points 
(7r/2, 0) and (7r/2, tt). Therefore the cyclotron mass scales 
as mo = V^/2, and the two sets of Landau levels are 
given by 

47r 1 

en{^) = t2(1 - r) ± —(n +-)(/; + •• • (13) 
When r = 1, the two sets of Landau levels separate (2.c). 

IV. SPECTRUM NEAR = 1/2 

We now concentrate on the region of the butterfiy spec- 
trum near cp = 1/2. Its peculiar structure can be ex- 




FIG. 5: Low field part of the Hofstadter spectrum e{(p), in 
the presence of a uniaxial potential, with r = 0.5; in the en- 
ergy range [—1,1]. The intermixing of the Landau levels is 
described in the text, and is given analytically by Eq.(13). 



plained from the spectrum precisely at (p = 1/2 which 
plays the role of a zero flux spectrum perturbed by a 
small magnetic effective flux / = — 1/2|, correspond- 
ing to an effective magnetic field Bf = 27r/ (the real 
field being n -\- Bj in our units). We first describe the 
low energy part of the spectrum. Then we study the 
vicinity of the band center where the spectrum consists 
in two Dirac cones which progressively merge when the 
staggered potential increases, until the value r = 1. 



A. Spectrum for cp = 1/2 

For cp = p/q = 1/2, the Hamiltonian has the form 

n{k) = - 

with r = As/2. The energy spectrum given by 

e{k) = ±2./cos2 ^ {cos ky-r)'^ (15) 



2{cosky-r) 1 + e , . . 

l + g2i/c. 2(r-cosA:^) ' ^ ^ 



is plotted on Fig. (6) for different values of r. It exhibits a 
pair of Dirac points and their merging occurs when r = 1. 



B. Spectrum near cp = 1/2, low energy 

An expansion of the dispersion relation (15) near its 
extrema at km = (0, tt) gives {q = k — km) 



e{km + q) = T2v/2 + 2r + r^±^lp^-t^ + . ■ ■ (16) 

V 2 + 2r + 





FIG. 6: Evolution of the electronic spectrum at p = 1/2^ 
when a uniaxial staggered potential (— l)^As is added. The 
four plots correspond respectively to r = As/2 = 0, 0.7, 1 and 
1.2. Note that that along the kx direction, we have represented 
two Brillouin zones, the first zone being [— 7r/2, 7r/2] . 



The cyclotron mass mo is therefore 



mo 



V2 + 2r + r2 



2^/l 



(17) 



so that in low "flux" / : 



(18) 



C. Spectrum near p — 1/2^ center of the band : 
Dirac points and their merging 

This study is the central point of our paper. When 
r = 0, the energy vanishes at the two inequivalent points 
= (7r/2,^7r/2) where ^ = ±1 is the valley index. In 
the vicinity of these two points, writing k = + g, the 
Hamiltonian can be expanded as 



niD^ + q)=2 



^qy -iQx 



(19) 



In each valley, the dispersion relation around is 

e(£>5 + q) = 2^q , (20) 

as can be seen on Fig. (6. a). 

We now apply a finite uniaxial staggered potential 
A5 ^ 0. The Dirac points move and are now located 
at the positions = (7r/2, ^ arccosr) (see Fig.6.b) . A 
second order expansion around these two points leads to 



n{D^^q) = 2 



-^Vl — r'^Qy — r — 



(21) 
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with a dispersion relation at small q of the form 



(22) 



with Cx = 2 and Cy = 2\/l — 

When r is increased, the two cones become anisotropic 
because the velocity along the y direction relating the two 
cones is reduced (Fig.6.b). When r = 1, the two cones 
have merged into a single one at the point = (7r/2, 0) 
(Fig.G.c), and the Hamiltonian becomes 



% 

2 



(23) 



leading to the "hybrid" dispersion relation 



e(fe) = ±2^^2+^4/4 . 



(24) 



Quite remarkably, this dispersion relation is linear in 
one direction, and quadratic in the other. The study of 
this semi-Dirac point has been the subject of a series of 
recent works^^^'^^ and the present problem constitutes 
a new realization of this merging. 

We now wish to describe the physics of the merging, 
using the effective Hamiltonian 



n = 



^ 2m* 

ic±q± 



-ic±q± 

^ 2m* 



(25) 



with q\\ = Qy and q± = q^- Within the rotation 
T-Lu =^1-L 7^"^ where 71 = 1/^2(1 - iay), the Hamil- 
tonian T-L is identical to the one Hu (see Eq.l) that 
has been recently introduced to describe the merging of 
Dirac points in a 2D crystal with time-reversal and in- 
version symmetries. When m*A > 0, there is a gap in 
the spectrum. Comparison between (25) and the expan- 
sion of (14) near Dq implies that A = 2(r — 1). When 
m*A < 0, the spectrum has two Dirac points separated 
by 2^—2771* A and the velocity cy near these Dirac points 

is C|| = 2A/m*. 

In order to quantitatively describe the merging of the 
Dirac points from the Hofstadter spectrum (r = 0) to 
the topological transition (r = 1), we choose to fix the 
velocities cy = Cy = 2\/l — and c± = Cx = 2. Then, 
the mass m* has to be fixed as m* = — 2A/cy = 1/(1 + 
r). The position of the Dirac points is given by q^ = 
=b2y^ slightly different from the real position q^ = 

zbarccosr (see ref. 9, for a discussion on this choice). 

Using the universal Hamiltonian (25) whose properties 
are known in a magnetic field, we now describe the evo- 
lution of the spectrum when approaching the merging. 

When r is small, the spectrum near e = can still be 
described as two independent cones with modified veloc- 
ities c_L = 2 and cm = 2\/l — r^- It is known that in 



this case,^*" the spectrum is quantized in two- fold degen- 
erate (due to the valley degeneracy) Landau levels with 
the dispersion relation = =b -y/2ncjcle5. In our nota- 
tions, this gives 



,(/)=±4(l-r2)i/4^/;^ 



(26) 



as confirmed on Figs.(7.a,b). When r increases, the do- 
main of validity of this expression is reduced. As seen 
in Figs. (2. b) and (7.b), the two- fold degeneracy of the 
Landau levels is removed when / is increased (in partic- 
ular the n = level), until, for r = 1, the degeneracy is 
completely removed, with a new field dependence of the 
levels. 

When r = 0, the Landau levels are doubly degenerate, 
due to the valley degeneracy of the two Dirac points. 
When r increases, the two valleys become coupled and 
the degeneracy is progressively lifted. The spectrum in 
the vicinity of (^^ = 1/2, can be described using semiclassi- 
cal arguments. It can be obtained from Bohr- Sommerf eld 
quantization 5 = 27r(n + 7)e5, where S is the area of a 
cyclotron orbit of energy e in reciprocal space. 

We have to distinguish two regimes (figure 6.b). In the 
low energy regime e < —A = 2(1 — r), this is the area 
enclosed by one of the two degenerate iso-energy lines 
encircling one Dirac point. For the Hamiltonian (1), this 
area has been calculated in ref. 9. Moreover, it has been 
argued that due to a Berry phase zLtt around each Dirac 
point, the mismatch factor 7 is 0. Therefore we find 



VT^{{e + A)K [r (J)] - AE [r (J)] } 



-nf 



(27) 



where K{x) and E{x) are complete elliptic integrals of 
the first and second kinds, respectively,^^ and R{x) = 
\/2x/ (x — 1). In the limit e <C —A, one recovers (26). 

In the high energy regime e > — A, that is above the 
saddle point, the area enclosed by an iso-energy line en- 
circles the two Dirac points so that, due to the cancella- 
tion of Berry phases, the mismatch factor 7 is now 1/2. 
The calculation of the area S{e) is this regime gives the 
semiclassical quantization rule 



{e + A)K 



R{e/A)\ 



2AE 



[R{e/A) 



]} 



= 37r2^/^T7(n' + -)/ 



(28) 



In the limit A ^ 0, one recovers (31). It is shown on 
Figs. (7), that this semiclassical approximation fits very 
well the numerical calculations, except in the vicinity of 
the line e = — A corresponding to the saddle point sepa- 
rating the two Dirac points (Fig. 6). 

To go beyond this semiclassical picture, we may explic- 
itly diagonalize the Hamiltonian (25) in a magnetic field. 
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Using the Landau gauge A\\ = Bx, and performing the 
Peierls substitution q\\ —eBx, one finds that the eigenval- 
ues Cn are solutions of an effective Schrodinger equation 
(see details in refs7'^) 



2/3 



{P^ ^ {5 ^ X^y - 2X)iIj (29) 



where P = {^^'q^ and X = (^)^ x 
are respectively dimensionless momentum and position 
([X, P] = i). We have introduced the dimensionless pa- 

/ * 2 2 \ 1/3 

rameter S = A/ f ^ ^^^^ j . In our case, this is a 

unique fonction of the flux / and of the parameter r 
describing the merging : 



.1/3 



(r-1) 



7^2/3(1 ^^)l/3y2/3 



(30) 



Eq. (29) is an effective Schrodinger equation for a par- 
ticle in a double weh potential V{X) = ((5 + X^)^ - 2X. 
It has been extensively studied in ref. 9, where its eigen- 
values are plotted as a function of the unique parame- 
ter S. Let us recall here its main characteristics. When 
^ << 0, the two wells of the potential are well separated, 
one recovers a \/rif spectrum of degenerate levels, prop- 
erly described by eq. (26), with an excellent fit of the 
butterfly spectrum near cp = 1/2, e = (Fig. 7. a). When 
—S is decreased, that is when r or / is increasd, the po- 
tential barrier between the two wells is reduced, and the 
tunneling between valleys becomes important. Therefore 
the two-fold degeneracy of the levels is removed, as seen 
on Figs.(7.b,c). 

When J = 0, that is r = 1, the spectrum exhibits a 
single "semi-Dirac" point. The effective potential V{X) 
is now a quart ic potential V{X) = X^ — 2X. In ref. 
3, a simple WKB quantization argument leads to e = 
±A{mc^Y^^[{n + l/2)cjc]^^^, where uOc = eB/m* and 
A = 3^/^7r/r(l/4)^/^. Using the values of the parameters 
mentioned above, we expect, at the critical point 
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e = ±G[{n'^l)f/' (31) 

with G = 2[367rVr(l/4)'^]i/3 7.99. This approxima- 
tion is quite good as soon as n > 1. However, it has been 
shown in ref. 3, that for the ground state, the prefactor 
has to be multiplied by a factor go :^ 0.808.^ Fig. (7. a) 
shows a remarkable agreement with the semiclassical cal- 
culation. 



V. CONCLUSION 



FIG. 7: Flux dependence of the energy levels of the tight bind- 
ing model on the square lattice, near cp = 1/2^ and for an 
alternate staggered potential characterized by r = As/2, a) 
The Hofstadter spectrum (r = 0). b) r = .5, and c) r = .75; 
the two-fold degeneracy of the levels is progressively removed, 
d) r — 1, at the merging of the Dirac cones, the Landau levels 
vary as [{n + 1/2)/^^'^. The continuous black lines are the 
results of the semiclassical quantization, which is quite good 
except, in the vicinity of the energy e = —A = 2(1 — r). At 
the merging r = 1, the energy levels are very well fitted by 
semiclassical calculation except the ground state for which a 
numerical factor has been introduced (see text). 



It is of great interest to study the physics of Dirac 
points, their motion and possibly their merging in con- lattice, the motion and merging has been known to be 
densed matter models. In the case of the honeycomb driven by a modification of hopping parameters, while it 
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is known that a modulated on-site potential opens a gap 
in the spectrum (like in the case of Boron Nitride), due 
to parity breaking. Here, we have shown a situation, 
the Hofstadter problem on a square lattice with half- 
flux quantum per plaquette (/)o/2, where the motion and 
merging of the Dirac points is not due to a change of the 
hopping integrals, but results from the application of an 
on-site uniaxial staggered potential. The merging of the 
Dirac points is discussed within a general hamiltonian, 
and the structure of the Landau levels near 00 /2, is ex- 



plained quantitatively. There has been a quite recent ex- 
citement about the possible realization of the Hofstadter 
spectrum with cold atoms in optical potentials.^^'^"*^ The 
merging of Dirac points described in this paper neces- 
sitates the application of a uniaxial staggered potential, 
which is possible with appropriate laser fields. 
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